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ABSTRACT 

We present a method for extracting the angular diameter distances, D^, and the expansion rates, 
H ^ of the universe from the two-dimensional Baryon Acoustic OsciUations (BAO) in the galaxy power 
spectrum. Our method builds upon the existing algorithm called the "fit-and-extract" (FITEX) 
method, which allows one to extract only D\/H from a spherically averaged one-dimensional power 
spectrum. We develop the FITEX-2d method, an extension of the FITEX method, to include the two- 
dimensional information, which allows us to extract Da and H simultaneously. We test the FITEX-2d 
method using the Millennium Simulation as well as simplified Monte Carlo simulations with a bigger 
volume. The BAOs, however, contain only a limited amount of information. We show that the full 
modeling, including the overall shape of the power spectrum, yields much better determinations of Da 
and i7, hence the dark energy equation of state parameters such as wo and Wa, than the BAO-only 
analysis by more than a factor of two, provided that non-linear effects are under control. 
Subject headings: cosmology : theory — large-scale structure of universe 



1. INTRODUCTION 

Dark energy, discovered via the observed luminosity 
dista i ices out to high-z Type la supernovae (jRiess et al.l 
Il998t iPerlmutter et al.lll999D. is the most mysterio us el- 
ement in physics today (see ICopeland et al.l [20061 for a 
recent review). 

As dark energy primarily affects the expansion rate of 
the universe, one can gain information on the nature of 
dark energy by measuring the cosmological distances as 
well as the expansion rates of the universe accurately.^ 

While the cosmic microwave background (CMB) and 
the Type la supernovae can be used for measuring the 
angular diameter distance out to z ~ 1090 and the lu- 
minosity distances out to z < 2, respectively, the power 
spectrum of matter distribution in the universe can be 
used to measure the angular diameter distances as well 
as the expansion rates of the universe out to a wider 
range of redshifts. 

Two length scales are encoded in the matter power 
spectrum, P{k) (see, e.g.. IWeinberd[2008l ): 

• The comoving Hubble horizon size at the matter- 
radiation equality, rnizeq) = c/[a{zeq)H{zeq)]. 

• The comoving sound horizon size at the so-called 
drag epoch at which baryons were released from 

photons, rs{zdrag) = /p^''''™''-' Cs{t)/a{t), where 
Csit) = c/ [^/3(l + a(i)3f7fc/(4f^^))] is the sound 
speed of photon-baryon fluid. 

The former determines the overall shape of the power 
spectrum of dark matter including the location of the 
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^ While dark energy also affects the growth rate of the am- 
plitude of matter fluctuations, which has been seen in the 
data via the so -called Integrated Sachs- Wolfe (ISW ) effect 
(e.g.,[Boughn fc C rittenden 2004; Nolta et al. 2004; Afshor dreFaLl 
120041) . we do not discuss the effect on the amplitude of fluctuations 
in this paper. 



peak of P{k) at k^q = l/r^f(zeg), whereas the latter de- 
termines the location of the baryonic features called the 
Baryon Acoustic Oscillations (BAOs). 

These length scales can be predicted from the 5- 
year data of the Wilkinson Micro w ave Anisotropy Probe 
(WMAP) ( Hinshawelall \200§. iDunklev et all [20081 : 
iKomatsu et al\\206^ : 
1 



^eq 



, • = (0.968 ± 0.046) X 10"^ Mpc" 

rH\Zeq) 

rs(2d™g) = 153.3 ±2.0 Mpc, 



(1) 
(2) 



and 
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ceq - ^-L,v,_;^5g, Zdrag = 1020.5 ± 1.6. (3) 

These lengths can be used as the "standard rulers," 
which give us the angular diameter distances as well 
as the expansion rates of the un i verse ([Seo fc Eise nsteinI 
[2001 iBlake fc Glazebrookl[200l iHu fc Ha imau 2003[).3 

We, as observers who measure the angular and redshift 
distribution of galaxies, can measure ]our distance ratios 
given by 



Bs{z)-- 



rnize 



il + z)DA{z) 
{1 + z)Da{z)' 



keqil + z)DAiz)' 



(4) 
(5) 



These predictions assume a flat universe and dark energy being 
the vacuum energy. For a non-flat universe with dark energy having 
a constant equation of state, w, the WMAP 5-year data yield k^q = 
(0.9751°:°^^) X 10-2 Mpc-i, rs{zarag) = 153.4^^;^ Mpc, z^g = 
1019.8 ± 1.5. 
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and Zjjr 



^ The matter power spectrum also contains the third distance 
scale, the Silk damping scale, which can also be used as the stan- 
dard ruler. The Silk damping scale is the smallest of these three 
distance scales, and its effect (i.e., the suppression of power below 
the Silk damping scale) is not as prominent as the effects of the 
other two distance scales. Nevertheless, the Silk damping must 
be taken into account when we model the full shape of the power 
spectrum. 
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rH{Zeq)H{z) _ H{z) 



keqC 



6Zs{z) : 



rs{zdrag)H{z) 



(6) 
(7) 



where Da{z) is the proper (i.e., not comoving) angu- 
lar diameter distance. We measure 9eq{z) and 0s{z) by 
comparing the predicted lengths with the corresponding 
observed lengths perpendicular to the line of sight, and 
5zeq{z) and 5zg{z) from the lengths parallel to the line 
of sight. ^ 

The BAOs have been detected in the current galaxy 
redshift survey data from the Sloan Digital Sky Survey 
(SDSS) and the Two-degree Field Galaxy Redshif t Sur- 
vey (MFGRS) (lEisensteiri et all [2005t iCole et all l2005l : 
lHiitsill2006t iPercival et al.ll2007D . However, the current 
data are not yet sensitive enough to y ield Da{z) and 
H{z) separately (lOkumura et al.l [20071 ): thus, one can 
only determine a combined distance scale ratio from the 
spherically averaged power spectrum. Since two spatial 
dimensions are available on the sky and one dimension 
is available along the line of sight, one can measure 



[el{z)5z,{z)] 



1/3 



[{l + zYD\{z)c/H{z)Y/^' 



(8) 



lEisenstein et al.l ([2005) have measured this quantity at 
z = 0.35 from the SDSS L uminous Red Galaxies (LRG), 
and IPercival et al.l ()2007[ ) have extended their analysis 
to include more data from the SDSS LRG, as well as the 
SDSS main galaxy samples and the 2dFGRS galaxies at 
z = 0.2. 

iKomatsu et al.l (|2008f) have combined these measure- 
ments with the CMB distance ratios determined from 
the WMAP 5-year data, the "WMAP distance priors," 
to obtain the constrai nts on dark ener g y prop erties. The 
analysis performed in IKomatsu et al.l (|2008f ) is a proto- 
type of what one can do in the future. It is clear that we 
can gain more information if we can measure Da{z) and 
H{z) simultaneously at various redshifts. Therefore, in 
the future we should be able to perform a much more sen- 
sitive test of dark energy properties by combining Da{z) 
and H{z) from the future galaxy survey data, and the 
CMB distance priors from the future CMB experiments 
such as Planck. 

Moreover, the BAOs capture only a part of information 
encoded in the shape of P{k). One would miss another 
baryonic feature, the Silk damping scale, by only measur- 
ing BAOs. A more serious drawback is that one would 
miss the other prominent standard ruler, fceg, completely, 
by only measuring BAOs. 

Nevertheless, there is one major advantage of us- 
ing BAOs: the phases (not the amplitude) of BAOs 
are less sensitive to the distortion of the shape 
of P{k) due to non-linear matter clustering, non- 
linear g alaxy bias, or non-lin e ar redshift space dis- 
tprtion (Seo & Eisenstein '2005'; 'Eisenstein et al.' '2007*; 
Nishimichi ct al. 2^007; Smith et al. 2008; Angulo et al. 
20081: ISanchez et af.ir2008l: ISeo et al.ll2008l ). As a resuh. 



* The measured power spectrum in redshift space is a function of 
the wavenumber parallel to the line of sight, fcy, and that perpen- 
dicular to the line of sight, fcj^, i.e., P = P(fc||,fcx). The angular 
observables, 9^q and ds, are measured from kj_, while the line-of- 
sight observables, Szeq and 5zs, are measured from 



many studies have focused on developing various ways to 
extract the distance information from BAOs. 

Most of the previous work focused only on ex- 
tracting the BAOs from the spherically average d 
P(k) (which giv es D\/H) (e.g.. IPercival et all [200% . 
lYamamoto et al.l (|2005i' ) have studied the monopole and 
quadrupole moments in the galaxy power spectrum and 
their implications for determinations of the dark energy 
equation of state parameter, w, and concluded that even 
in the worst case scenario (i.e., absence of the BAOs fea- 
ture on the observed power spectrum), galaxy survey can 
still provide useful limits on w from a combination of 
the monopole and quad. r upole power spectra. Recently, 
iPadmanabhan fc White! (|2008f ) have explored an extrac- 
tion of the quadrupole moment of the two-dimensional 
power spectrum, P(fc, //), which gives a different distance 
combination, DaH. 

In this paper, we shall develop a method for extract- 
ing Da and H simultaneously from the two-dimensional 
BAOs. Since we do not use spherical averaging or trun- 
cate the Legendre expansion of BAOs at arbitrary or- 
ders, our method uses more information than most of 
the previous methods. To our knowledge, the full two- 
dimensional extr action of Da and H fr om BAOs has been 
explored only bvlWag ner et al.l (|2007f ). 

This paper is organized as follows. In §|2]we give a brief 
account of the original one-dimensional "fi t-and-extract" 
FITEX) method, which was developed bv lKoehler et al.l 
2007t l for extracting BAOs from a spherically averaged 
one-dimensional P{k). We then extend this method to 
the two-dimensional FITEX-2d method by including the 
full two-dimensional information without spherical aver- 
aging. In §|3|we extract Da and H from simulated noisy 
data using the FITEX-2d method, and show that the 
FITEX-2d yields unbiased estimates of Da and H. In 
§|4|we repeat the same analysis for a more realistic simu- 
lation , using the Millennium Simulation (ISpringel et al.l 
I2005D . In §|5]we propagate errors in H{z) and DAiz) 
to those in the dark energy equation of state with the 
parametrization of w{z) = wq + Waz/^l + z). We con- 
clude in § ini 

Throughout this paper we shall use the cosmologi- 
cal parameters given by Q,n — 0.277, f^A = 0.723, 
Qb = 0.0459, Hs = 0.962, an d h = 0.702 punklev et all 
120081 IKomatsu et"al] l2008l ). which are the maximum 
likelihood value s inferred from the WMAP 5-year data 
([Hinshaw et al. ' 2008) combined with the current BAO 
data (IPercival et al. 2007,,) and Type la supernova data 
(jKowalski et al.ll2008D ~ 

2. FITEX-2D: METHODOLOGY 

We develop a method for extracting Da and H si- 
multaneously from the two-dimensional BAOs without 
spherical averaging. 

Our method builds upon the exi sting "fit-and-e x tract" 
(FITEX) method developed by iKoehler et all (pOOTi ) 
for extracting D\/H from a spherically averaged, one- 
dimensional P{k). The FITEX method extracts BAOs 
by fitting and removing the non-oscillatory part of P{k), 
which leaves only th e oscillatory component, i.e., BAOs. 
IKoehler et all (|2007D model the non-oscillatory, smooth 
part by the following functional form: 



P. 



Id 

smooth 



{k) = 



A 



1 + Bk^ 



(9) 
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where Us is the primordial tilt, while A, B, 



and a are free parameters. iKoehler et al.l (|2007l ) 
have shown that this function is flexible enough to 
fit out the smooth part of the spherically averaged 
P(k) measured fro m the Hubble Volume Simulation 
(jEvrard et all l2002l ). They have tested the FITEX 
method particularly for a large scale, k < 0.3 h Mpc~^, 
at high redshifts, 1.9 < z < 3.8, that arc relevant to 

the Hobby Eberly Dark Energ y Experiment (HETDEJC; 

Hill G. J. and Gebhardt. K. and Komatsu. E. and MacQueen. P. J. I 

mm . 



where /3 = f/bi is a fc-independent function that de- 
pends on the linear galaxy bias, bi, and the cosmological 
parameters (mainly Qm) via 



dlnD 
dlna ' 



(16) 



where D is the growth factor of linear density fluctua- 
tions. We therefore find 



We make a simple extension of the one-dimensional 
FITEX method by including angular dependence. We 
model the two-dimensional smooth power spectrum by 



l + g^^\k)P^{^i) 



smooth 

(k) 

g^'\k)p,i^,) 



9 



(6) 



{k)Pe{fi)]lO) 



where is the cosine of the angle 9 between k and the line 
of sight, i.e., fi = cos 9 and tan 9 = k±/k^^. Therefore, 
jj, = and /i = 1 for fcy = and k± = 0, respectively. 
Here, Pi{i-t) is the Legendre polynomials: 



P2(p)=\(:3f^^-i), 

-P4(/i) = i (35^*^-30/^2 + 3), 



(11) 
(12) 



^6(m) = 



1 

16 



231/ -315/ + 105^2 -5). (13) 



The odd multipoles must vanish by symmetry. One may 
include Z > 8 if necessary, but we find it sufficient to 
include the terms only up to ^ = 6. 

Finally, ^('^(fc) is given by the 6th-order polynomials 
with only even powers of k: 



g('\k) 



(14) 



where all of af's are varied simultaneously for each I. 
The odd p owers must van ish because they are not ana- 
lytic in k (| Weinberg! l2008l ) . We include the terms only 
up to fc^, as we include the multipoles up to Z = 6. If, for 
instance, Z = 8 is included, then fc* may also be included 
for consistency. 

Aside from the primordial tilt, rig, the FITEX-2d con- 
tains 17 free parameters (5 for Pj^oothi^) P^^^ 4x3= 12 
for the angular dependence). While it may sound like 
many, the number of data points available on the two- 
dimensional power spectrum is usually much larger, and 
thus our fit is well behaved. 

It may be instructive to use the conventional model for 
the redshift space power spectrum to show what these 
parameters are supposed to capture. The leading order 
angular distortion is given by the so-called Kaiser effect, 
which arises from coherent convergi ng velocity f low to- 
ward the linear overdensity region (iKaisedHosl . The 
linear Kaiser power spectrum is given by 



PkaTserik, m) = &f (1 + 2/3/ + /3>^)P 



■ylinear 



(fc) 



.hi 



2 1 

1 + 77/3 + 7/3 



pi' 



(15) 



1/3(1 + 1/3) 



l + i/3+i/32 



,(4). 
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l + |/3+i/32 



(17) 
(18) 

(19) 



and the other terms are zero. 

Another example is the so-called Finger-of-God (FoG) 
effect, which arises from random motion within virialized 
halos. When the distribution of the pairwise peculiar 
velocity within a halo is given by an ex ponential dis- 
tribution with the vel ocity disper sion a'}. (IPeeblej 119761 : 



iDavis fc Peebleslll983D . one finds (|Ballinger et al.lU996f ) 

:(A:,/ 



PFoG{k,^l) 



kaiser ' 

1 + pB^i-^al ■ 



(20) 



While the FoG yields many terms when expanded into 
the Legendre polynomials, it is still a good approxima- 
tion to truncate the expansion at Z = 6 if fc is sufficiently 
smaller than l/tr^. Note that the FoG effect yields terms 
in the form of powers of (fc/i)^; thus, it makes sense to 
use the same number for the maximum power of fc (see 
Eq. (fT4|) ) and the maximum multipole (see Eq. (fTO| ) of 
the FITEX-2d fitting function. 

In general, neither of these two expressions are ade- 
quate. The linear Kaiser formula is valid only on very 
large scales, while the exponential FoG formula is valid 
only on very small scales. At the intermediate scales 
we find more complicated ex pressions from, e.g. , the 
3rd-order perturbation theory (jHeavens et al.l[l998h . To 
account for these complications we have included fc- 
dependent coefficients for the Legendre polynomials. 

In Figure [1] and [2] we show the performance of 
P^mooth^k^ l^)- Figure [1] we show a simple analyti- 
cal model^ for the non-linear galaxy power spectrum in 
redshift space given by 



Fg(fc,/ = 6? [F55(fc) + 2/3/F5(,(fc) +/3V^'ee(fc)] 
1 



1 + /2fc2/i2^2 ' 



(21) 



where Pss{k), Psg{k), and Peeik) are the density- 
density, density-velocity, and velocity-velocity power 
spectra computed from the 3rd-order perturbation the- 
or y, and they are gi ven by Eq. (63), (64), and (65) 
in IScoccimarrol (l2004f). respectiyely. This form is simi- 
lar to Eq. (71) ofl coccimarr 3 (l200l . but we have re- 
placed exp(— /2fc2/i2(j2^ g^jjj j- }jjg formula by 1/(1 -|- 
/2fc2/t2j5.2-j respectively, where = 0.6cr2 is the 

1-d peculiar velocity dispersion with an empirical fudge 

^ This model is admittedly too simple to be realistic. We shall 
test the FITEX-2d method in a more realistic setting using the 
Millennium Simulation in §[4] 
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Fig. 1.— Illustration of the FITEX-2d method. This figure shows 
an anisotropic non-linear galaxy power spectrum before we apply 
FITEX-2d. The contours show ln[P(A:|| , fcx)] at z = 2, where we 
have computed P(A:||,A:j^) from Eq. I I21I I. Anisotropic distribution 
of power due to redshift space distortion is apparent. 

factor of 0.6 calibrated o ff our simulations presented in 
iJeong fc Komatsul (|2006f ). Here, is given by 

We chose z = 2 and bi = 2.5. The contour of power 
spectrum is anisotropic in Fig. [1] due to the redshift space 
distortion; however, we recover isotropy after subtracting 
the best-fitting Psmoothi^ ^ f^) from the anisotropic data 
(see Fig. [2]) . We see that the BAOs have been extracted 
successfully, with isotropy of the oscillation phases recov- 
ered well. 

3. EXTRACTION OF Da AND H FROM NOISY DATA: 
FITEX-2D VS FULL MODELING 

In § 13.11 we show how well we can estimate Da and 
H from the two-dimensional BAOs extracted from noisy 
data using the FITEX-2d method. In 13.21 we compare 
the BAO results to the accuracy one would obtain from 
the full modeling of P(fc, /^), including the overall shape. 
In other words, for the former (BAOs) we only use 9^ and 
6zs for measuring Da and H, while for the latter (full 
modeling) we can use Og, dz^, Oeq, Sz^q, as well as the 
Silk damping scale for measuring Da and provided 
that non-linear effects (non-linear matter clustering, non- 
linear redshift space distortion, and non-linear bias) are 
under control. 

Note that the treatment of non-linear effects in this 
section is too simple to be realistic. For a more realistic 
treatment we shall use the gala xy power spectru m from 
the Millennium Simulation (Spr ingel et alJl2005f ) in § H 

3.1. FITEX-2d 

To estimate errors in Da and H from the FITEX-2d 
method, we use simple Monte Carlo simulations. 

For the underlying spectrum we use the same data as 
shown in Fig. [1] which includes a simplified modeling 
of non-linear matter clustering and non-linear redshift 
space distortion as given by Eq. (HI]). As for the galaxy 
bias, we use a linear bias with bi ~ 2.5. 

Once the underlying spectrum is specified, it is 
straightforward to compute the errors in Pg{k\^,k±), 




[h Mpc"'] 



Fig. 2.— Illustration of the FITEX-2d method. This figure shows 
the power spectrum shown in Fig. [T] minus the best-fitting two- 
dimensional smooth spectrum, P^^^joth^^W ' ^-l)' given by Eq. mUl l. 
The structure of BAOs, i.e., the oscillatory feature, is now appar- 
ent. The FITEX-2d method recovers the isotropic distribution of 
the BAO phases successfully, which makes it possible to use the 
distribution of the phases for measuring and H simultane- 
ously. (Top) Positive BAO peaks. (Bottom) Negative BAO peaks 
(troughs). 

(TPg, provided that the distribution of Pg{k^\,kj_) is 
a Gaussian. We use the standard formula that in- 
cludes sampling varianc e as well as shot noise (see, e.g., 
IJeong fc Komatsull2008[ ) 

Q'Pfllfclhfcx) _ / 1 l + ngPg{k\\,k±) 

Pg{k\\,k±_) \j Vsurveyk±Ak±Ak\\ ngPg{k\\,k±_) 

(23) 

where Ug is the number density of galaxies, Vsurvey is 
the survey volume, lS.k±_ and Afcy are the fundamental 
wavenumbers, i.e., the resolution in fcj^ and fc||. We take 

1 /3 

these to be Afc|| = A/c^ = {2it) /Vsurvey 

We use fjpg from Eq. ((23| to calculate the r.m.s. er- 
ror in Pg{k\\,k±), and generate 1000 Monte Carlo re- 
alizations. We then apply the FITEX-2d method to re- 
move the smooth component from each realization to ex- 
tract BAOs. For each realization, we measure Da and 
H simultaneously by fitting the phases of extracted two 
dimensional BAOs to those of the reference BAOs ex- 
tracted from either (i) the linear power spectrum, or (ii) 
the non-linear power spectrum given by Eq. (|2ip . with 
known DA,ref and Href- (Later we find that using the 
linear spectrum as the reference BAO yields the biased 
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estimates of DA,ref and Href-) We use a simplex down- 
hill method for ^^-minimization in the two-dimensional 
parameter space. The number of free parameters for this 
analysis is two, i.e., Da and i7, and we do not include 
the amplitude in the fit. We have checked that including 
the amplitude does not change the results very much, as 
the amplitude and the phases of BAOs are nearly un- 
correlated (see Appendix IA.3.11 for more details) . This 
is true in both real and redshift space. When we apply 
FITEX-2d to the simulated data, we perform a fit out to 
= 0.40 /i Mpc-^ 



We choose the survey parameters, Vs 



survey j 



and 



1.10 



1.00 



0.95 



1.05 



1.00 



such that they roughly match those expected for the 

Hobby-Eberly Dark Energy Experi ment (H ETDEX) 

(iHill. G. J. and Gebhardt. K. and Komatsu. E. and MacQueen. P. Jo.fe s 
12Q04D: Ng 



with 



0.755 

coverage 

-3 ri„„3 6 



10^ 
of 



and 1.9 

420 deg2, 



< z < 3.5 
which yields 



the sky 

^survey — 3.0 h GpC 

We iind that, when the phases extracted by FITEX-2d 
are compared with the reference BAOs extracted from 
the linear power spectrum, the best-fitting values of Da 
and H averaged over 1000 simulations disagree with the 
underlying, "true" values by 0.05% and 0.63% for Da 
and H, respectively, due to the phase shift of BAOs 
caused by non-linearities (including non-linear redshift 
sp ace distortion) . This res ult extends the previous study 
bv lNishimichi et al] (j2007f ). who studied a spherically av- 
eraged 1-d power spectrum and found that the bias was 
less than 1% in {D\H~'^ f/^. 

On the other hand, when the phases are compared with 
the reference BAOs extracted from the non-linear power 
spectrum (Eq. (HT])), the best-fitting values of Da and 
H agree with the true values to well within the Monte 
Carlo sampling error; thus, we confirm that the FITEX- 
2d method yields unbiased estimates of Da and H . 

In Figure [3] we show the projected error ellipses on 
Da and H from the BAOs extracted with the FITEX- 
2d (larger, dotted contours; same in all four panels). 
We find 1.8% and 2.5% errors on Da and respec- 
tively, with the cross-correlation coefficient of r = 0.44, 
from the Monte Carlo simulations. For the same sur- 
vey parameters, the B AO Fisher matrix proposed by 
iSeo k Eisens"teinl (l200l yields 1.5% and 2.5% errors on 
Da and H, respectively, with r ~ 0.41. Therefore, we 
conclude that the FITEX-2d method yields the results 
that nearly saturate the Fisher matrix bound, i.e., it is 
nearly an optimal method in a sense that it can yield 
the smallest errorbars one can obtain with the BAO-only 
analysis. 



3.2. Full Modeling 
To calculate the errors in Da and H expected from 



the full modeling of the two-dimensional galaxy power 
spectrum, Pa(k,ii), we use the Fisher matrix given by 
(see, e.g.. lEisenstein et al.|[l99l ISeo fc Eisensteinll2003l) 



(27r)3 



^ ^^d\nPg{k,fi) d\nPg{k,^i) 



0.90 
0.90 



InA 


InA, /5 

■ ■ 

^-Full modeling 

BAO fitting 


InA, /S, a^^ 

■ :@ 


InA, /S, CT ^, n , a 



0.95 1.00 
D /D 
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D /D 



Fig. 3. — Accuracy of Da and H extracted from BAOs with the 
FITEX-2d method appHed to simulated Monte Carlo realizations 
that appro xima te the HETDEX survey (the larger, dotted con- 
tours; see ^ I3.1|l . The best-fitting values of Da and H agree with 
the true values; thus, the FITEX-2d method yields unbiased esti- 
mates of Da and H. The solid contours show Da and H from the 
full modeling, including the overall shape of the power spectrum, 
with various parameters marginalized over. (Note that the BAO- 
only contours are unaffected by the marginalization.) For this we 
have used the Fisher matrix forecast (see § 13.21 1. The inner and 
outer ellipses show 68% and 95% C.L., respectively. (Top Left) 
the full modeling Fisher matrix is marginalized over the overall 
amplitude, InA, (Top Right) marginalized over InA and the linear 
redshift distortion parameter, 13, (Bottom Left) marginalized over 
InA, /3, and the velocity dispersion in the FoG factor, ct^, (Bottom 
Right) marginalized over InA, /3, cr^, and the shape of the initial 
power spectrum, and as- 



where 9i = (In Z?yi, IniJ, In A, /?, tj^, tIs, as) for i = 1, 
2,.. .,7, respectively, kmax = 0-4^0 h Mpc^^, where A is the 
overall amplitude of the power spectrum, /3 is the linear 
redshift distortion parameter, ct^ is the calibrated 1-d 
velocity dispersion (see Eq. ((2T|) ). and Ug and as describe 
the shape of the initial (primordial) power spectrum: 



(25) 



Here, the weight function, u;(fc,/i), is one half of the 
so-called "effective volume," 



w{k,y) 



1 + ngPg{k,fi) 



-Veff{k,^l). 



(26) 

The effective volume is equal to the actual survey volume, 
^survey: in the Sampling variance dominated regime, 
Pg{k,fi) ^ whereas it is small in the shot-noise 



09,^ 



^ The HETDEX is expected to detect 0.755 million Lyman-a 
emitting galaxies between 1.9 < z < 3.5 over 420 deg^ in 3 years 
of observations on the Hobby-Eberly Telescope. 



dominated regime, Pgik,^) <C ^/rig. In Figure |4] we 
show UgPgi^k, ^) for Ng — 0.755 x 10^ and bi — 2.5 as a 
function of z. The factor of 1/2 accounts for symmetry 
in k ^ — k. The derivatives of lnPg(/c,/i) with respect 
to 9i are calculated and given in the Appendix IA.2I 
Unlike for BAOs, which are insensitive to the parame- 
w{k, /i), ters that affect the overall shape, for the full modehng we 
need to make sure that we take into account potential de- 
generacy between Da and H and any other parameters 
that affect the overall shape. In this paper we include 
\nA, (3, ci\,^, rig, and as- (We shall comment on the 
effects of non- linear bias in S 13.31) . 



(24) 
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TABLE 1 

Fisher matrix forecast for, the full power spectrum analysis with various choices of 

marginalization 



none In A (3 a"^ ris cts In A, /3 

A In Da (%) 0.279 0.877 0.317 0.282 0.479 0.416 1.100 

AlnH {%) 0.437 0.786 1.124 0.801 0.509 0.539 1.134 

tda h 0.382 -0.720 -0.309 0.082 -0.227 -0.226 0.038 

Alnfl'(%) 0.187 0.762 0.317 0.259 0.386 0.363 0.775 

P, crl In A, al In A, (3, cr'^ In A, Us, as In A, jS, Us In A, f3, a^, 

Os Wa, as 

A In Da (%) 027 0891 LlOl L089 033 050 

AlniT(%) 1.457 1.101 1.468 0.984 1.362 1.530 

H -0.383 -0.632 0.005 -0.820 -0.199 -0.098 

Alnii'(%) 0.322 0.869 0.879 0.974 1.000 1.014 

Note. — The fractional errors in Da and H, and their cross-correlation coefficients, 
rDA,H, and the fractional errors in the combined 1-d distance scale, R (Eq. (|29|) ). marginal- 
ized over several combinations of para meters: In^, /3, g^, Us and Us . The cosmological pa- 
rameters are taken from Table 1 of the lKomatsu et al.l (120081 ) ( "WMAP-I-BAO-I-SN ML"). 
The survey parameters approximate those of HETDEX: the survey area and target redshift 
are 420 deg^ and 1.9 < z < 3.5, respectively, the number of galaxies is Ng — 0.755 x 10^, 
and the bias is assumed to be linear with bi = 2.5. 
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Fig. 4. — The galaxy power spectrum times the number density 
of galaxies, ngPg{k, fj,), where the number of the galaxies is fixed 
for each re dsh ift bin to Ng = 0.755 X 10^, and Pg{k, ^) is computed 
from Eq. I I21I I with bi = 2.5. The shot noise dominates the error 
budget when ngPg{k, fi) < 1. Contour values are [0.1, 0.3, 0.5, 1.0, 
3.0, 5.0]. {Top Left) z = 1, {Top Right) z = 2, {Bottom Left) z = 3, 
{Bottom Right) z = 4. 

We study the effects of marginalization over various 
parameter combinations by taking the submatrix, Fij, 
of the full 7x7 matrix with the index, i, oi di running 
from 1 to 7, such that the submatrix includes the matrix 
components of desired parameters to be marginalized. In 
other words, the parameters that are not included in the 
submatrix are fixed and not marginalized over. 

Then, we compute the marginalized errors in In Da 
and In H as 

ai„D. = V'(i^-i)ii, (27) 

CTi„ff = y^(F-i)22. (28) 

To simplify the analysis, we fix all the other cosmo- 
logical parameters, such as /(z), fib/i^, etc. These cos- 
mological parameters will be determined by the future 



CMB mission, Planck, accurately, and therefore it is a 
good approximation to simply fix them, and vary only 
XiiDa, \nH, A, f3, a^, Ug and as- The fiducial value for 
the bias is set to bi = 2.5 and / — d\nD/d\na is com- 
puted from the fiducial cosmological model.'' We expect 
that the analysis of the bispectrum (Fourier transform 
of three-point function) will give a precise determination 
of bi (as well as non-hnear bias parameters such as 62) 
(jSefusatti fc Komatsul l2007f ). and therefore it is also a 
good approximation to simply fix it. However, we also 
explore a more conservative case where we do not know 
what bi is, i.e., we marginalize over the overall amplitude 
as well as /?. In the future work we also plan to inves- 
tigate the effect of marginalization over 62, using a joint 
analysis of the power spectrum and bispectrum. There- 
fore, our calculation presented here will provide the lower 
limit to the errors in In Da and IniJ expected from the 
full modeling of the power spectrum measured in a sur- 
vey like HETDEX. We use the same survey parameters 
that we have used in § 13.11 and we integrate Eq. ([24]) up 

to krnax = 0.40 h Mpc"^ 

In Figure [3] we show the resulting error ellipses from 
the full modeling, in the smaller, solid contours, with 
four choices of marginalization. (We present the results 
from more choices of marginalization in Table [T]) First, 
for all choices of marginalization we find that the sizes 
of the errors in both Da and H are substantially smaller 
than those from the BAO-only analysis with the FITEX- 
2d. For example, determinations of both Da and H are 
improved by more than a factor of two in the case of 
the amplitude marginalization. This is expected, as we 
are able to use more information encoded in the power 
spectrum; namely, the Hubble horizon at the matter- 
radiation equality epoch and the Silk damping scale. 
Second, Da and H are anti- correlated for the amplitude 
marginalization, with the cross-correlation coefRcient of 

^ One might also wish to marginalize over / for the following 
reason: while / can be calculated from the cosmological parameters 
assuming the validity of General Relativity, one might choose to 
let / free and use it for testing the validity of General Relativity. 
In this paper we chose to assume the validity of General Relativity, 
but one can extend our analysis to let / free in a straightforward 
manner. 
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r — —0.72 (see top-left panel of Fig. [3|), as opposed 
to a positive correlation seen in the BAG-only analy- 
sis. This is due to the marginalization over the overall 
amplitude: if we fixed the overall normalization, then we 
would still find a positive correlation between Da and H 
with r = 0.38. 

The origin of the negative correlation is the so-calle d 
Alcock-Paczynski (AP) test (jAlcock fc Paczvnskiiri979[ ): 
when the redshift space distortion is known perfectly 
well, the departure of the power spectrum in redshift 
space from isotropy, i.e., dependence of P{k,^) on /i^, 
can be used to determine DaH, resulting in r = —1 for 
a power-law power spectrum. The contributions from 
departures of P{k) from a pure power-law, i.e., the ex- 
istence of "standard rulers," such as BAGs, the Hub- 
ble horizon at the matter-radiation equality and the 
Silk damping scale, make r bigger than —1. (See Ap- 
pendix IA.3I for more details.) When In^ and f3 are 
marginalized over simultaneously, the correlation be- 
tween Da and H nearly disappears: the AP test no 
longer works when we marginalize over the linear red- 
shift space distortion. We find r = 0.038 (see top-right 
panel of Fig. [3|) for this case. 

When In A is marginalized over while the other pa- 
rameters (/?, (7^, Us, and Us) are held fixed, we find 
0.88% and 0.79% errors on Da and H, respectively, with 
r = —0.72. The more parameters we marginalize over, 
the greater the cross-correlation coefficient between Da 
and H as well as the errors on Da and H become. Note 
that the increase in the errors does not necessarily imply 
the decrease in the statistical power in constraining dark 
energy properties: since the cross-correlation coefficient 
is also reduced, the error in the combined 1-d distance 
scale, R, is much less affected by the marginalization 
(see Table [T]) . The error in In R has been computed as 
(jSeo fc Eisenstein .2007): 



'In Da 



(1 



1 + 2rain Da /^In H + < /crf„ ' 



(29) 



Finally, the errors in Da, H, and R for various choices 
of marginalization are: {<Jin d a t '^in h , o'ln r) = (0.88%, 
0.79%, 0.76%), (1.10%, 1.13%, 0.78%), (1.10%, 1.47%, 
0.88%), and (1.25%, 1.53%, 1.01%) for the marginaliza- 
tion over In^, In^ and (3, In A, (3 and a^, and In^, /3, al, 
ris and as, respectively (see Table[T]for more comprehen- 
sive list). This result should be compared with that from 
the BAG-only analysis: (crin - c^in if, crin i?)=(1.76%, 
2.47%, 1.08%). It is clear that the full analysis, even 
with a generous set of marginalization choices, beats the 
BAG-only analysis with a significant gain in the distance 
determination accuracies. 

3.3. Caveat for the full modeling 



Gur analysis presented in § 13.21 is too simplistic and 
optimistic, as it ignores any systematic errors due to 
our lack of understanding of the effects of various non- 
linearities in the power spectrum. 

Among the three major non-linearities, non-linear mat- 
ter clustering is under control, at least for high redshifts, 
i.e., z > 2, as one can model non-linear evolution of 
matter fiuctuation s almost exactly by the 3 rd-order per- 
turbation theory (jJeong fc Komatsul 120061 ). While the 
nominal 3rd-order perturbation theory breaks down at 



lower redshifts, z ~ 1, there have been a number of stud- 
ies aiming at improving upon our ability to compute P(k) 
at z ^ 1 or even low er redshifts ( Crocce fc Sco ccimarr 
l2008t iMatarrese fc P ictroni 200 71: iTaruva fc H iramats' 



3^ 



2008'; l Valageasll200'7[ IMatsubarall2008l : lMcDonaldl l2007i 
Therefore, it is quite possible that the non-linear matter 
clustering will be fully under control in the near future, 
at least for the scales that are relevant to the BAG scales, 

i.e., k < 0.40 h Mpc~^ 

In a separate paper (jJeong fc Komatsull2008f ). we show 
that non- linear galaxy biasing is also under control in the 
weakly non-linear regime. Gne can use the perturbation 
theory approach combined with th e local bias assumption 
(|Frv fc Gaztanagal 119931 : iMcDon ald 2006) to model the 
galaxy power spectrum with non-linear bias. 

The most problematic one is the non-linear redshift 
space distortion. Gur understanding of non-linear red- 
shift space distor tion, especially the Finger-of-God (FoG) 
effect, is limited (|Scoccimarrci||200^ . Therefore, whether 
one can achieve the accuracy of Da and H {H in particu- 
lar) reported in Fig. [3] depends crucially on our ability to 
correct for the FoG effect. This is work in progress. Note 
that the marginalization over tr^ should capture some of 
the increase in the errors in distance scales due to our 
ignorance of FoG. 

4. EXTRACTION OF Da AND H FROM THE MILLENNIUM 
SIMULATION 

How realistic is our result for the determinations of 
Da and H from the BAG phases using the FITEX-2d 
method? Since our Monte Carlo simulations used in 
§ 13.11 arc too simple, in this section we test the FITEX- 
2d method further b y using the Millennium Simulation 
(|Springel et al.ll200l . 

We use the Millennium Galaxy catalogue, generated by 
the s e mi-analytical galaxy formation code (iBower et al.l 
I200I IBenson et alll200a iCole et al. ' 2000). We have 
measured the two-dimensional power spectrum of galax- 
ies in redshift space from the Millennium Simulation, and 
applied the FITEX-2d method to remove the smooth 
component. We then find the best-fitting Da and H from 
the BAG phases extracted from the FITEX-2d. Again, 
we use the data up to kmax = 0.40 h Mpc~^. 

In Figure [5] we show the result. The best-fitting 
value that we find from the Millennium Simulation cor- 
responds to one point at the center of the contours. 
We find the errors from the Monte Carlo simulations 
that we described in § 13.11 with the survey parame- 
ters replaced by those of the Millennium Simulation: 
Vsurvey = (0.5 /i"! Gpc)^, Ug = 0.138 Mpc'^ , and 
z = 3.06. (There are 17,238,935 galaxies in the Millen- 
nium Simulation at z — 3.06.) For the theoretical power 
spectrum that we use for generating Monte Carlo sim- 
ulations, we use the best-fitting power spectrum for the 
galaxy catalogue of the Millennium Simulation found in 
Jeong fc Komatsu (2008). 

Since the volume of the Millennium Simulation is ~ 24 
times as small as that would be surveyed by HETDEX, 
the uncertainties in Da and H are larger for the Mil- 
lennium Simulation. (Compare Fig. [5] with the larger 
contours of Fig.[3l) We find 5.1% and 6.8% errors on Da 
and H , respectively, with the cross-correlation coefficient 
of r = 0.43, from the Monte Carlo simulations. These er- 
rors are larger than those from HETDEX Monte-Carlo 



8 



SHOJI, JEONG & KOMATSU 




Fig. 5. — Accuracy of Da and H extracted from BAOs with the 
FITEX-2d method appUed to the Millennium Galaxy Simulation 
in redshift space a t z = 3 ((S pringcl et al. 2005; Bower ct al. 200^: 
IBenson et al.l2003l : [C"ole et a l. 2000). The best-fitting values oi Da 
and H agree with the true values to within statistical errors of the 
Millennium Simulation; thus, the FfTEX-2d method also yields 
unbiased estimates of Da and H for the Millennium Simulation. 
The solid and dotted lines show 68% and 95% C.L., respectively. 



simulation by a factor of two (rather than \/24 ^ 5) as 
the shot noise on the power spectrum of the MiUennium 
Simulation is much smaller than that of HETDEX sim- 
ulation. 

The best-fitting values of and H are well within 
68% C.L. region, which indicates that the FITEX-2d is 
able to yield unbiased estimates of the BAO phases from 
the Millennium Simulation. 

These results indicate that the FITEX-2d method that 
we have developed in this paper can be used for extract- 
ing the BAOs and measuring Da and H safely from the 
real data. It would be interesting to apply the FITEX- 
2d method to the two-dimensional powe r spectrum mea- 
sured from the SDSS LRG samples (jOkumura et al.1 
|2007() . and extract Da and H from them. 

5. ERROR PROPAGATION TO THE DARK ENERGY 
EQUATION OF STATE 

In § [3] and § |4l we have estimated errors in Da and 
H from two different approaches, i.e., the BAO fitting 
using the FITEX-2d method and the full modeling. In 
this section, we propagate errors in Da and H to those 
in the dark energy equation of state parameters. We 
parametrize w(z) using the linear model, w(z) = wn + 
Waz/{l + z) (lLindeill2003HChevallier fc Polar skill 2 OOll) . 

We propagate the errors in Da and H to those in wo 
and Wa by 



P _ dp, dp J ^ 



(30) 



where is the Fisher matrix for the dark energy 

parameters, Fij is the Fisher matrix for Da and 
Pi = (In 13^, In iJ) for i — \ and 2, and = (wo,u'a) 
for CL — \ and 2. 
Partial derivatives of Da and H with respect to wq 
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Fig. 6. — Partial derivatives of In Da and ln_ff with respect 
to the dark energy equation of state parameters, uiq and uia, 
as a function of z for two different cosmological models. (Left) 
(wo, Wa) = (-1.0,0.0). (Right) {wo,Wa) = (-1.1,0.5). 

and Wa are given by 
din Da 



z')f{z')g{z')-^/^dz' 



dwo 
dh-iDA 

dWa 



3o ilMl 

2 ^ J^^ g{z')~^f^dz' 



= — Ha 
2 



J^^[\n{l + z')~j^]f{z')gizr'/'dz' 



d\nH 3 



J,' g{z')-'/'dz' 



5W 



d\nH 3 



dWa 



Ha 



ln(l + z) 



1 + z 



/(f), 



(31) 

,(32) 

(33) 
(34) 

(35) 
(36) 

Figure [6] shows the derivatives as a function of z be- 
tween 0.5 < z < 6.5 in two different cosmological mod- 
els, (wo, Wa)=(-l-0, 0.0) and (-1.1, 0.5). The for- 
mer is the ACDM model, while the latter resembles 
the maximum likelihood v alues of wg and Wg fr om the 
WMAP-hBAO-HSN-hBBN (jKomatsu et al.|[200l) . The 
derivatives are similar for these cosmological models, and 
therefore we use the ACDM model as the fiducial model 
for computing the derivatives. 

We add the distance information from CMB as 



where f{z) and g{z) are given by 



/(z)=exp 



5(2;) = ri„ 




htotal/\ f^CMB 
^al3 (Z) = i'af) 



(37) 



where we assume that the CMB experiment yields 1% 
determination of the angular diameter distance out to 
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Fig. 7. — Projected 68% constraints on the dark energy param- 
eters, wq and Wa- the BAO fitting with the FITEX-2d method 
(dotted) versus the full modeling (solid). For both cases, we use 
the power spectrum up to kmax = 0.40 h Mpc~^, and we assume 
that the CMB experiment measures the angular diameter distance 
out to 2 = 1090 with 1% accuracy. The survey area and the num- 
ber of galaxies are 420 deg^ and Ng = 0.755 X 10® for all cases. 
(Top Left) 0.5 < 2 < 1.5, (Top Right) 1.5 < z < 2.5, (Bottom Left) 
2.5 < z < 3.5, (Bottom Right) 3.5 < z < 4.5. 



1090, 



I.e., we use 



1090) d\nDA{z ^ 1090) 



dqa 



(38) 

We are interested in how the BAO-only analysis com- 
pares with the full modeling. In Fig. [7] we show the pro- 
jected error contours on wq and Wa calculated from the 
BAO-only analysis with the FITEX-2d and those from 
the full analysis at four redshift bins: 0.5 < z < 1.5, 
1.5 < z < 2.5, 2.5 < z < 3.5, and 3.5 < z < 4.5. The 
survey area and the number of galaxies are 420 deg2 and 
Ng = 2.9 X 10^ for all redshift bins. From the BAO-only 
analysis we find (Aw^Awa) = (0.29,1.26), (0.38,1.39), 
(0.55, 1.92), and (0.91,3.18), whereas from the full mod- 
eling we find {Awo,Awa) = (0.09,0.27), (0.06,0.17), 
(0.09,0.35), and (0.17,0.68), for 0.5 < z < 1.5, 1.5 < 
z < 2.5, 2.5 < z < 3.5, and 3.5 < z < 4.5, respectively. 

We therefore conclude that the full analysis yields 
much better constraints on wq and Wa than the BAO- 
only analysis. 

6. CONCLUSION 

In this paper we have developed a method, called 
the FITEX-2d method, to extract the two-dimensional 
phases of BAOs from galaxy power spectra in redshift 
space. Our model builds on and extends the existing 
one-dimensional algo rithm, called FITEX, developed by 
IKoehler et all ((2007h . 



Our method removes the smooth, non-oscillating com- 
ponent from the observed galaxy power spectrum in 
redshift space. The fitting function consists of the 
smooth one-dimensional spectrum that depends only on 
Plioothi^) given by Eq. multiplied by the angle- 
dependent function expanded in the Legendre polynomi- 
als with even multipoles. The coefficients of the Legen- 
dre polynomials contain even powers of k. The resulting 
function, given by Eq. (fTO|) . is able to capture the non- 
oscillating part of the galaxy power spectrum well. 

We have tested the FITEX-2d method using the ana- 
lytical model without any noise, the Monte Carlo realiza- 

tions with noise expected from the HETDEX experiment 

(iHill. G. J. and Gebhardt. K. and Komatsu. E. and MacQueen. P7 
l2004f ). and the galaxv catalogue created from the Mil- 
lennium Simulation (Springel et al. 2005). In all cases 
the FITEX-2d method yields unbiased estimates of the 
angular diameter distance. Da, and the expansion rate, 
H. 

However, the BAOs capture only a part of distance 
information encoded in the galaxy power spectrum. To 
exploit the distance information, especially the equal- 
ity scale, rH{zeq), we have explored the constraints on 
Da and H from the full modeling of the galaxy power 
spectrum in redshift space. Provided that three key 
non-linearities (non-linear matter clustering, non-linear 
galaxy bias, and non-linear redshift space distortion) are 
under control, we find that the full modeling yields the 
constraints that are better than the BAO-only analysis 
by more than a factor of two both in Da and H , and the 
dark energy parameters such as wq and Wa- 

While the effects of non-linear matter clustering 
(iJeong k K omatsu' 120061; ICrocce & Scoccimarrd l2008|; 
i Matarrese fc Pictroni 2007; Taruva & Hiramatsu 200i; 
IValageasI |2007i; iMatsubara 2008: McDonald 2 007) and 
non-linear galaxy bias ([jcong & Komatsu 200^ are be- 
ing understood in the weakly non-linear regime that is 
relevant to the future galaxy surveys, the effects of non- 
linear redshift space distortion are poorly understood. 
While the FITEX-2d method that we have developed in 
this paper are useful for obtaining robust constraints on 
Da and H, hence the dark energy properties, one must 
understand non-linear redshift space distortion to fully 
exploit the full information content of the galaxy power 
spectrum in redshift space. We would then be able to 
reduce the errors in Da and H by more than a factor of 
two. 



We thank the anonymous referee for illuminating com- 
ments on the correlation coefficients, which motivated 
our doing more thorough analysis. This material is based 
in part upon work supported by the Texas Advanced Re- 
search Program under Grant No. 003658-0005-2006. EK 
acknowledges support from an Alfred P. Sloan Research 
Fellowship. The Millennium Simulation databases used 
in this paper and the web application providing online 
access to them were constructed as part of the activities 
of the German Astrophysical Virtual Observatory. 



10 SHOJI, JEONG & KOMATSU 

APPENDIX 
FISHER MATRIX CODE 

In this Appendix we describe what we have implemented in our Fisher matrix code, which is pubUcly available 
as a part of "Cosmology Routine Library (CRL)," developed by one of the authors (EK). This code includes the 
non-linear matter power spectrum in both real and redshift space, as well as marginalization over the amplitude, the 
linear redshift space distortion, the velocity dispersion of Fingers-of-God, the primordial tilt and running index. In 
the future release we plan to include non-linear galaxy bias and primordial non-Gaussianity. 

Basics 

A simple, approximate formula of the Fisher matrix for galaxy survey is given by fe.g.. lSeo fc EisensteinI [20031) 



AnPdk d\nPg{k,fi) d\nPg{k,fi) 

'WW Jo ^ 

where Pg{k,^) is the galaxy survey power spectrum calculated theoretically as a function of parameters, 6i are the 
parameters to be extracted from the data, and /i) is a function given by 



ngPg{k,fi) 



^survey i,-^^') 



_l + ngPg{k,^) 

Here, Ug and Vsurvey are the number density of galaxies and the volume of survey, respectively. 
In linear theory, Pg{k,fi) is given by 

Pgik,fi)=blR{fi')P'^"'^^^k), (A3) 

where bi is the scale independent linear bias factor, p'""^'"'^/^) jg the linear matter power spectrum, and R{iJ.^) describes 
the linear redshift space distortion effect (Kaiser effect): 

i?(/i2)^(l+/3;,2)2 (A4) 

/3=(dln£)/dlna)/5i, (A5) 
where D is the growth factor of the linear density fluctuations, and a is the scale factor. 

Derivatives 

To calculate the logarithmic derivatives of P{k) in Eq. (|A1[) . let us write down the non- linear galaxy power spectrum 
(with linear bias) as (Eq. ([21]) ): 

Pg{k,ii) = hi [Pss{k) + 2(]^i^Pseik)+(3^fi''Pgeik)] x ^ ^ y2fc2^2^2 ■ 

We compute the derivatives with respect to the following seven parameters: the angular diameter distance. Da, the 
Hubble expansion rate, H, the overall amplitude of the galaxy power spectrum, A, the linear redshift space distortion 
factor, f3 = f/bi, the velocity dispersion with an empirically calibrated fudge factor, ct^, the tilt of the primordial 
power spectrum, n^, and the running index, as {Pini oc in[fe/fcpi„ot]-j^ choose the convention such that 

) = {\nDAMH,\nA,p,dl,ns,as). (A6) 



'1, C2, (73, (74, (75, (76, (77 
- 2 



The derivatives with respect to In A, /?, (t^ , , and as are easy to evaluate. They are given by 

91nFg(fc,M) 



1, (A7) 

(A8) 



d\nA 

dlnPg{k,fi) _ 2fi^Pse{k) + 2|3^i*Pee{k) 



813 Pss{k) + 2Pfi^Pse{k) + p'^^i^Pee(k) 

d\nPg{k,^i) _ 

ddl ""l + /2A:V^2' 

din Pg{k,^) _ dlnPiniik) 



(A9) 



dug dus 
dlnPgik,fi) _d\nP,^,{k) _ 1 



Ink, (AID) 



In 



kpivot 



das das 2 

We compute the derivatives with respect to InD^i and ln_ff in a two step process. First, we write 

d\nPg{k,n) _ d\nPg{k,^i) d\nk dlnPg{k,fi) V 

dlnDA dink din Da ^ df? cflnD^' 

dlnPg{k,^i) _ din Pg{k,^i) dink dlnPg{k,^i) 9/^2 
dlnH ^ dink dlnH 'd^T^ dlnH' 



(All) 

(A12) 
(A13) 



where 



ainfc 

din Da 
dink 
dlnH 

_dil_ 

5 In Da 
9/^2 
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(A14) 
(A15) 

-2a*'(1-A*'), (A17) 



dlnH 

dlnPg{k,y) 2l3Pse{k) + 2l3''ti'Pee{k) pP^l 



d^i^ Pss{k) + 2pfi^Pse{k)+f3^tJ-^Pee{k) ^ + pk^^l^al 

Finally, we need to know the "effective spectral index", Uef f{k, ii), given by 



(Af8) 



dlnPo(k. n) , . 

ne//(fc,/i)^ ainV ' ^^^^^ 



or explicitly 



neff{k,^i) = 



Pss{k)nss{k) + 2l3^i'Pse{k)nse{k) + p'^ ^l^Pee{k)nee{k) 
Pssik) + 2Pfi^Pse{k) + f3^fi^Pee{k) 

2pe^?di 



(A20) 



where 



(A22, 

"'^'^(^^^ ginfc ■ ^^^^^ 
Correlation Coefficients 

In this subsection we explore the behaviour of the cross-correlation coefficient between and H in various cases. In 
particular we focus on the effect of the marginalization over the overall amplitude with f? IA.3.2p and without f? IA.3.ip 
the additional marginalization over the redshift space distortion. 

No redshift space distortion, /3 = 

Let us evaluate the Fisher matrices in the limit that the redshift space distortion is absent, i.e., /3 = 0. In this 
limit, the weighting function in Eq. (jAl|l and the effective spectral index in Eq. (|A19p become independent of fi, i.e., 
w{k,^) w{k) and Ueffik^fj,) Ueffik). We obtain 

^11 - -^[n^fjmMk) - 1^''? (A24) 

k'^dk , 



Fi2=/ -^[n,ffik)]Mk) I d^^il-^^^){-^,') (A25) 



Fi3- -^n,ff{k)w{k) dfi{l-f,') (A26) 

F22^ -^[nefs{k)Yw{k) d^,{-^?f (A27) 

^ne//(fc)t.(fc) d^,{-^,^) (A28) 

F^, = j^ ^^ik)J^dt, (A29) 

Now, in order to understand the effect of the structure of n^ff, let us assume that the galaxy power spectrum is a 
pure power-law, i.e., neff{k) = n and n is the independent of k. In this limit, we obtain 
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8n^ 


2n^ 


2n 


15 2 


i5 

n 


3 




n 




5 


3 




n 


1 


3 


3 



F^,=^\ ^ I (A30) 

where w = J ^^^w{k). 

The marginahzed errors of parameters and the correlation coefficients are computed from the inverse of the Fisher 
matrix. However, one can show that the matrix given in Eq. (jA30p is singular. In other words, Da and H are 
completely degenerate with the amplitude for a power-law power spectrum. This result shows that only the departure 
of the power spectrum from a pure power-law, i.e., the existence of characteristic scales, can break the degeneracy 
between Da and H, and A. These scales are often called the "standard rulers." 

To understand the structure of the Fisher matrix in Eq. (|A30[) better, let us add small perturbations, e > 0, to the 
diagonal elements, and invert the matrix. The result is 

1 1 n 

I (2+n2)e (2+ri2)e (2+n^}e 

{F-% = \ -j2lWe (2+^ (2^^ \+0{e") (A31) 

n n 

^ {2+n'^)e (2+n^)e (2+n^)e 



We find that the correlation coefficient between Da and H is 



^ ■■ - -1 (A32) 



V(F-i)n(F-i) 



22 

as £ — *■ 0. Therefore, In Da and InH are totally anti-correlated, which implies that, although we cannot determine 
ln_Dyi and IniJ simultaneously, we can determine In Da -I- IniJ = In(Z)ytJf), even for a power-law power spectrum. ^ 
This is known as the Alcock-Paczyhski (AP) test (jAlcock &: Paczvnskilll979[ ). 

There is a special case in which the covariance between A and Da or H may be ignored. One may imagine the 
situation where neff{k) depends upon k such that A is uncorrelated with Da or H. For example, if neff{k) oscillates 
about zero, then J^"^'"' liA^rief f{k)w(k) would be small compared with J^''""" ^{j^y^^yjij^') or 

Therefore, Fi^ and ^23 may be ignored, making A de-correlated with Da and H. In this case, the Fisher matrix is a 
2-by-2 matrix given by 

Fzj = wn^ ( A ~F ) (A33) 

V 15 5 / 

The inverse of this matrix is then 

1 / 9 3 \ 

{F-% = — i I (A34) 
The correlation coefficient between Da and H is thus given by 

ri2 = , ^/^ = ^ - 0.408. (A35) 
v/9/4 X 6 VQ 

This result has been derived bv lSeo &: Eisensteinl (|2007[ ). and justifies the use of BAOs as a way to measure Da and 
H with a correlation coefficient of 0.408. 
From these studies we are led to the following conclusion: 

• When the information is dominated by BAOs, the correlation coefficient between Da and H is ri2 ~ 0.408. The 
amplitude of the BAOs contributes little to the errors on Da and 7J, as the amplitude information is de-correlated 
with Da and H. 

• When the information is dominated by the AP test, ri2 ~ — 1. 

• In reality, as we have shown in this paper, BAOs contribute less than the overall shape of the power spectrum. 
Also, the shape of the power spectrum is not exactly a power-law. As a result, the correlation coefficient from 
the full analysis is usually negative (or small positive), but always greater than —1 (see Table [T]). 

With redshift space distortion, f3 > 

Next, let us consider the case where the redshift space distortion cannot be ignored. In this case, we see from 
Eq. (|A2p and Eq. (|A20p that the weighting function, w{k, n), and the effective spectral index, neff{k, fj), are no longer 
independent of /i. The analytical treatment is also possible for this case, although the results are too complicated to 
be useful. We therefore report on the numerical results. 

Here, we choose the survey parameters given in 13.21 with the non-linear power spectrum of Eq. (|2ip . The results 
from the numerical calculations of the Fisher matrix are given in Table [T] We find that the marginalization over the 

* The other cross-correlation coefficients are ri3 ^ =pl and r-^s ±1 for n > and n < respectively. 
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Fig. A1. — Projected 68% constraints on the dark energy parameters, wq and Wa. The full modeling (solid) marginalized over different 
combinations of parameters as well as the BAO-only analysis (dotted) are shown. For all cases, we use the power spectrum up to 
kmax = 0.40 h Mpc~^, and we assume that the CMB experiment measures the angular diameter distance out to z = 1090 with 1% 
accuracy. The survey area and the number of galaxies are 420 deg2 and Ng = 0.755 X 10^, and the redshift range is 1.9 < z < 3.5 for all 
cases. {Top Left) marginalized over InA, {Top Right) marginalized over InA, f3, [Bottom Left) marginalized over InA, /3 and tr^, {Bottom 
Right) marginalized over InA, /3, ct^, rts and as- 

amplitude information, and that over the amplitude and the shape of the primordial power spectrum (i.e., Hs and as) 
give the cross-correlation close to —1; thus, one relies on the AP test. The marginalization over the amplitude and 
the linear redshift space distortion (i.e., /3) drive the cross-correlation towards zero, as the AP test no longer works 
when the linear redshift space distortion is marginalized over. However, in both cases the errors in the combined 1-d 
distance scale, R, are about the same. In other words, while one changes the orientation of the ellipse, the area is 
approximately preserved. 

In summary, when the amplitude information is marginalized over, the information is mostly coming from the 
dependence of P(fc,/i) on which yields a constraint on DaH via the AP test, while when both the amplitude and 
the linear redshift space distortion are marginalized, the most information is coming from the standard rulers, which 
can constrain Da and H separately, driving the cross-correlation towards zero. 

Finally, in Fig. lAll we show how different choices of marginalization over parameters influence the error contours of 
Wo and Wa- {AwQ,Awa) = (0.08,0.27), (0.08,0.30), (0.24,0.85), and (0.24,0.86), for the cases of no marginalization, 
marginalization over InA, marginalization over InA, /3 and trj, and marginalization over InA, (3, ct^, Us and 
respectively. 

User's Guide 

When using the Fisher matrix code, one may choose the form of the model galaxy power spectrum from: 

• Linear power spectrum with the linear redshift space distortion (Kaiser effect), 

• Non-linear power spectrum from the 3rd-order perturbation theory with the linear redshift space distortion, 

• Non-linear power sp ectrum from t he 3r d-order perturbation theory with the non-linear redshift space distortion 
given by Eq. (71) of iScoccimarrol ()2004[ ). or 

• Non-linear power spectrum from the 3rd-order perturbation theory with the non-linear redshift space distortion 
given by Eq. (PTjl . 

Next, specify the number of parameters one wishes to marginalize over, and then choose the parameters from: InA, 
/3, crj, ns, and as. 

A given galaxy survey can be sliced up into multiple redshift bins. After entering the survey area in units of deg^, 
one is asked to enter the following parameters at each redshift bin: the redshift range {z„iin < z < Zmax) , the number 
of galaxies in units of millions in the bin, &i, k^ax in units of h Mpc~^, and the redshift error in units of km/s. 

Th e linear power spectrum at z = 30 has been precomputed using the CAMB code (jLewis, Challinor fc Lasenbvl 
I2000D for the maximum hkelihood parameters given in Table 1 of iKomatsu et all (|2008f) ("WMAP-hBAO-hSN"). The 



14 



SHOJI, JEONG & KOMATSU 



ingredients of the non-linear power spectra, Ps6^ Pse, and Pee, have been precomputed from the Hnear spectrum at 
z — 30. These spectra are then evolved to a specified redsh i ft by the appropriate growth factor obtained by solving 
the differential equation given in Eq. (76) of iKomatsu et all (|2008f ). 

Finally, the code yields the errors on In Da, \nH, rDA,H, and Ini? (see Eq. (p9)) for the definition of the error in the 
combined distance scale, K). 
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